Mon. Not. R. Astron. Soc. 000,[I]i5](2007) Printed 2 February 2008 (MN MgX style file v2.2) 



Bianchi Vllh models and the cold spot texture 

M. Bridges^*, J. D. McEwen\ M. Cruz^, M. P. Hobson\ A. N. Lasenby\ 
P. Vielva^, E. Martmez-Gonzalez^ 

^Astrophysics Group, Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CBS OHE, UK 

^Instituto de Fi'sica de Cantabria, CSIC-Universidad de Cantabria, Avda. los Castros s/n, 39005, Santander, Spain 



Accepted — . Received — ; in original form 2 February 2008 



ABSTRACT 

We have returned to our previous Bianchi Vllh analysis in light of the ICruz et all (l2007bh 
suggestion that the cold spot observed near th e southern Galactic pole may be a remnant 
temperature perturbation of a cosmic texture. In lBrid^es et alJ (l2007l) we found two favoured 
left handed Bianchi Vllh templates with restricted prior probabiUties so that the template was 
centred close to the cold spot. Using WMAP data 'corrected' for the texture fit we have now 
reexamined both models to assess any changes to these conclusions. We find that both mod- 
els are left almost entirely unconstrained by the data and consequently exhibit significantly 
reduced Bayesian evidences. Both models are now disfavoured by the data. This result rein- 
forces our previous suggestion that the cold spot was likely to be driving any Bianchi Vllh 
detection. 

Key words: cosmic microwave background - methods: numerical - methods: data analysis. 



1 INTRODUCTION 



The anomalous cold spot discovered dVielva et alj2004l : ICruz et aP 
I2OO5L I2OO6L l2007ah in Wilkinson Microwave Anisotropy Probe 
(WMA P ) observations of the cosmic microw ave background 
(C MB) teennett et al.ll2003l : iHinshaw et aDl2007h has been shown 
bv ICruz et al. I l l2007bl) to be consistent with a temperature pertur- 
bation induced by a cosmic texture. A cosmic texture is a particu- 
lar type of cosmic defect predicted by certain theories of high en- 
ergy physics. They form at symmetry breaking phase transitions in 
the early Universe and are extremely energetic events producing 
both hot and cold spots in the CMB. The texture hypothesis is the 
most convincing explanation of the cold spot anomaly made to date, 
however evidence for the texture requires further observational sup- 
port. One of the most promising aspects of the texture hypothesis is 
that tests may be performed using future observations of CMB po- 
larisation, which may either substantiate or refute evidence for the 
presence of a texture. The discovery of a texture would have such 
profound implications for our understanding of the early Universe 
that the cold spot-texture hypothesis certainly warrants further in- 
vestigation. 

Bianchi Vllh models in which the universe exhibits a global 
rotation and shear have also been considered in an attempt to ex- 
plain many of the anomalies in the WMAP data. A positive de- 
tection of a Bian chi Vllh ter nplate embedded in the WMAP data 
was first made bv I Jaffe et"al] (|2005). After 'correcting' the WMAP 
data for this template a number of previously reported anomalies 
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in the data disappear 1 Jaffe et al. I I2OO5I : iLand & Magueiid bood 
ICavon et al.l2006| - |McEwen et alj2006h . However, the correspond- 
ing best-fit Bianchi V\\ templa te was shown to be incompati- 
ble wi th concordance cosmology jjaffe et al.ll2006l) . [Bridges et al.l 
(l2007h performed a more rigorous Bayesian MCMC analysis to 
determine the evidence for Bianchi Vllh models, concluding that 
there is weak evidence against a Bianchi Vllh component when the 
axis of the Bianchi coordinate system is allowed to vary over the 
entire sky. However, weak evidence/or a Bianchi template remains 
when the axis is restricted to lie in the direction of the cold spot. 
These results suggest that it may have been the cold spot that wa s 
driving the positive template detection made bv lJafi^eetal] ( l2005h . 

The focus of this work is to determine whether there is evi- 
dence for any Bianchi Vllh component embedded in the WMAP 
data once the d ata is 'corrected' fo r the cold spot texture template 
determined by ICruz et alJ l l2007bh . If all positive evidence van- 
ishes, and one accepts the texture explanation of the cold spot, then 
Bianchi Vllh models may be rejected definitively. The remainder 
of this letter is structu red as follows. Firstly, the procedure followed 
bv lCruz etaLH2007bl) to fit a texture template to the cold spot in the 
combined, foreground-cleaned Q-V-W map (hereafter WMAP co- 
added map) is discussed, before the template is used to 'correct' 
the WMAP internal linear combination (ILC) map (since the sub- 
sequent MCMC analysis requires a full-sky map and the co-added 
map necessitates a Galactic cut). Using different processed versions 
of the WMAP data for template fi tting and analysis is acceptable 
and was shown bv iMcEwen et alj (|2006) not to alter analysis re- 
sults for the Bianchi case. Secondly, the Bayesian evidence for the 
texture 'corrected' WMAP ILC map is computed and compared 
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to the evidence computed previously dBridges et alj|2007h for the 
original data. Concluding remarks are then made. 



2 COLD SPOT TEXTURE FITTING 



L oc e T", 



The cold spot has been shown bv lCruz et al.l j2007bb to be consis- 
tent with a temperature perturbation induced in the CMB by cos- 
mic texture at redshift z ~ 6. Unwinding events associated with 
texture induce cold and hot spots in the CMB. An analytic approx- 
im ation for the temperatur e profile produced by a texture is given 
bv dXurok & Spergelll 19901) 



(1) 



^ - +g 

where 6 denotes angular separation, 9c is the scale parameter of 
the texture and e is the amplitude parameter related to the symme- 
try breaking scale (f>Q. This temperature profile is an approximation 
and is not valid for large co-moving scales. Consequently, for prac- 
tical applications the profile is truncated beyond its half-maximum 
by matching its val ue and derivative to a Gaussian, as discussed 
bv lMagueiid l ll995l) . 9c is then equal to the standard deviation of 
the matching Gaussian. The resulting texture temperature profile 
was fitted to the cold spot by Cruz et al.l |2007b) and found to be 
favoured to the null hypothesis of no texture. I n this section we re - 
view the texture fitting procedure performed bv lCruz et alj ( l2007bh 
and present the best-fit texture profile determined. 

The co-added m ap of the three-year WMAP data release 
dHinshaw et alj l2007h is use d to fit the texture profile to the 
cold spot. Cruz et al. I ( l2007hl) d egrade the co-adde d map in the 



HEALPi)(J pixelization scheme jGorski et alj 1200^ to a resolu- 
tion parameter of A'sidc = 64. This resolution is sufficient to re- 
tain all information associated with the cold spot, which occurs 
on a half- width scale of ~ 5°, and reduces the number of pixels 
used in the template fitting, thereby reducing the computational 
cost of the fitting procedure. The template fitting was performed 
in a circular area of 20° radius centred at Galactic coordinates 
{b = -57°,/ = 209°). Although the total angular size of the cold 
spot is ~ 10°, it is necessary to consider a 20° radius patch in or- 
der to take into account the entire neighbourhood of the spot since 
the original Sph erical Mexican H at Wavelet analyses that high- 
lighted the spot jVielva et al.ll2004l) convolves all pixels in this re- 
gion. These pixels could therefore contribute in an important way to 
the detected structure and must be included when fitting the texture 
template. 

A hybrid Bayesian-freq uentist approach is considered for the 
template fitting performed bv lCruz et al ]| |2007bh . The Bayesian ev- 
idence ratio is used to perform hypothesis testing, which is then 
calibrated using Monte Carlo simulations. The data are found to 
favour the alternative hypothesis that a texture is present. Pa- 
rameters e and 0c of the temperature template profile T, defined 
by l[T}, are then determined by maximising the posterior probabil- 
ity for the alternative hypothesis, where the posterior is given by 
Bayes' Theorem: 



Pr{<d\D, H,) oc Pr(D\&, H,)Pr(&\H.,), 



(2) 



with likelihood Pr(D\&, and prior Pr{&\H,,), where D repre- 
sents the WMAP co-added data and = (6, 9c). The likelihood 
function is assumed to be Gaussian: 
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where 



(D-TfN"'(D-T), 



(3) 



(4) 



with the generalised noise covariance matrix N including both 
CMB and noise contributions. The calculation of the noise con- 
tribution to N is straightforward since the noise in the WMAP co- 
added map is uncorrelated and well defined through the number 
of observations per pixel. In order to obtain the CMB contribu- 
tion to the matrix, the covariance function for the WMAP co- 
added map is calculated taking into account both pixel and beam 
effects. As a complementary test the CMB cova riance matr i x is cal- 
culated through 70,000 Gaussian simulations. ICruz et sl\ ( l2007bl) 
compared the values computed from the data to those obtained 
using simulati ons and found error s were negligible. As a conserva- 
tive prior on 6, ICruz et all (l2007bh choose < e ^ 10 '', th e COBE- 
normalised amplitude jPen et alJl9"94l : lDurrer et aljl999h . A scale- 
invariant distribution of texture spots on the sky is predicted 



cWr 



3 9i 



(5) 



where y is a dimensionless constant and k a fraction of unity. In or- 
der to obtain the prior for the scale parameter 0c > IS is normalised 
to unity between t/mi,, and 0,nax- Photon diffusion would smear out 
textures smaller than a degree or so, hence the lower 9c bound is set 
to ^^^lin = 1°, according to the resolution considered. At large scales 
textures ar e rare, hence the u pper 9c bound is set to S^ax = 15°. In 
this setting Icruz et al ] l l2007bl) obtain texture parameter estimates 
of 6 = 7.7 X 10"' and 9c = 5.1°. The original co-added data, the 
best-fit texture template an d their difference ar e shown on a small 
patch on the sky in Fig 1. of lCruz et al ] (l2007bh . 



3 EVIDENCE FOR BIANCHI VII* MODELS 

Bianchi VII), models induce characteristic 'spiral' temperature fluc- 
tuations on the CMB sky. They are described by four independent 
quantities: a matter energy density n„,, a dark energy density Q.^, 
the current vorticity a) and h, which physically relates the char- 
acteristic wavelength over which the principle axes of shear and 
vorticity change orientation, and determines the 'tightness' of the 
spiral pattern. The pattern position is defined by Euler angles a, p 
and 7 13 Additionally one must specify the direction of rotation or 
han dedness of the spir al . 

Ijaffe et alj ( l2005h found evidence of a significant correlation 
between a Bianchi Vllh model and the WMAP 1-year data. In 
Bridges et al. (2 007h we performed a more rigorous Bayesian anal- 
ysis that confirmed the iJaffe et al. I (j2005) fit in WMAP 3-year 
data. In this framework our entire inference is contained in the 
multidimensional posterior distribution from which we can extract 
marginalised parameter constraints and the comparative Bayesian 
evidence to select the most appropriate model parameterisation. We 
aimed to use the evidence to establish whether it was necessary to 
include a Bianchi Vllh component in addition to a standard ACDM 
cosmology. We concluded that the only significant Bayesian evi- 
dence favouring such an inclusion existed where the central posi- 
tion of the Bianchi Vllh spiral was fixed close to the cold spot. In 

^ We adopt the active zyz Euler convention con'esponding to the rotation of 
a physical body in a fixed coordinate system about the z, y and z axes by y, 
/} and a respectively. 
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Table 1. Summary of Bianchi Vllh component priors used in this analysis. 



Full Bianchi 

nl, = [0.01,0.99] 
= [0.01,0.99] 
h = [0.01, 1] 
0) = [0,20] X 10"'" 
r = [0, 2;r]rads 
ChiraUty = L/R 



Table 2. Cosmological and Bianchi parameterisations for each of the pa- 
rameter subsets studied. 



Model 


Cosmology 


Bianchi 


A 


A, 




D 


A, 


nf„„ h, oj, r, L/R 


G 


As 


Q.^,h,oj, y, L/R 



this letter we aim to establish how the signific ance of thi s conclu - 
sion has changed in light of the texture fit of ICruz et al.l ( l2007bh . 
We use the WIVIAP 3-year IL C map (Fig. [T] (a)) subtracting the 
cold spot template described bv lCruz et al.l ( l2007lj) (Fig. [T](c)) to 
yield the 'corrected' map shown in Fig.[T](b). 

The models we will consider here are those of A, D & G from 
lBridgesetal1 ( l2007h . for convenience we summarise the parameter 
combinations and priors in Tables [T] & [21 IVIodel A is a standard 
ACDIVT cosmology, with one free parameter As, the amplitude of 
scalar fluctuations, this being the only parameter constrained on 
these large scales. Both models D and G had the Bianchi template 
axis fixed close to the cold spot via heavily constrained priors on 
a and p. However D and G differ in the inclusion of a dark en- 
ergy density component which significantly opens up the parame- 
ter space along a d egeneracy in the Qa - ^^m plane (see Fig. 6 in 
iBridges et al.l2007h . It was first suggested by Jaffe et al. (2006) that 
a dark energy term might give the Bianchi Vllh models more free- 
dom to accomodate a close to flat cosmology preferred by WIvlAP 
and others. Although the Bianchi degenera cy found is intriguing ly 
close to the geometric degeneracy seen in Spergel et al.l ( l2007h it 
does not cross the WMAP confidence contours above the 2cr level. 
From this we concluded that the set of cosmological parameters 
preferred by WMAP are incompatible with those of any Bianchi 
Vllh component. We proceeded in our analysis by decoupling the 
Bianchi Vllh matter and energy densities from their ACDIVI coun- 
terparts so that the Bianchi VIIi, component was essentially treated 
as a template with energy densities Q.^ and Cl^. Left-handed mod- 
els D and G were the only ones found to show marginally signifi- 
cant evidence favouring their inclusion and so we have chosen these 
to carry over for this analysis. IVIorphologically both D and G are 
very similar (owing to the degeneracy described above); the best 
fitting t emplate is shown in F ig.[T](d). 

In Bridges et alj ( l2007h both models D and G resulted in 
highly constrained marginalised posterior distributions (red lines of 
Fig.[2]c&[3} on each of the Bianchi Vllh parameters. In this analysis 
(black lines of Fig.lJl&O however all parameters except perhaps 
a slight preference for h ~ 0.2, are entirely unconstrained, and cru- 
cially, a non-zero likelihood is observed in the vorticity al a> = 0. 
In the Bianchi Vllh formalism a> is highly correlated with Bianchi 
Vllh signal amplitude so such a result illustrates that with the tex- 
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-0.300 ^^^^^^^^^^^^K .^^^^^^^^^^M +0.300 

(a) WMAP ILC map 




(b) Cold spot con-ected WMAP ILC map 




-0.300 ^^^^^^^^^^^^^^^^^V^^ ^^^^^^^^^^^^H +0.300 



(c) Cold spot template 




-5.000E-02 ^^^^^^^^^^^^^^^^^^H. ^^^^^^^^^^M .5.000E-0? 



(d) Previous best-fit Bianchi template 

Figure 1. Full-sky CMB and correction maps displayed in the Mollweide 
projection. The original WMAP ILC map is illustrated in panel (a), while 
the cold spot corrected map is shown in panel (b). The cold spot template 
is displayed on a full-sky ma p in panel (c). In pan el (d) the best-fit Bianchi 
Vllh template determined bv lBridges et alj (2003) is shown. Notice that the 
cold spot texture template aligns closely with the central cold swirl of the 
best-fit Bianchi template. All maps are given in units of mK. 
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Table 3. Bayesian evidence differences (logarithmic) between models A 
and left-handed D & G for the original WMAP 3-year ILC map and that 
connected for the texture fit of Cruz et al. ( 2()07b). 



Data \ Model 


D(L) 


3-year WMAP 


+ 1.2 ±0.2 


3-year WMAP (texture corrected) 


-0.1 ± 0.2 


Data \ Model 


G(L) 


3-year WMAP 


+ 1.2 ±0.2 


3-year WMAP (texture corrected) 


-0.8 ± 0.2 




0.2 0.4 g 0.6 0.8 



tot 

Figure 2. Marginalised Bianchi VIIi, parameters for left-handed model D 
usin g the original WMA P ILC map (red) and that corrected for the texture 
fit of lCruz et alj j2007bl) (black). 



ture corrected map the data no longer prefers any additional compo- 
nents of this kind over a standard cosmology. Although redundant, 
given the unconstrained posteriors, the same conclusions can be 
confirmed via the Bayesian evidence estimates (see Table[3j. Both 
models D and G record significantly reduced evidences, model D 
is now about as favoured as a ACDM cosmology alone (model A) 
while model G i s now, almost s ignificantly, disfavoured when com- 
pared with A. In lBridges et al.l (2007) all right-handed models were 
found to have disfavouring evidences. It is worth noting that this 
result remains unchanged with the texture corrected data. Both left 
and right handed models are now equally disfavoured. 




4 5 6 7 8 O-^ 0-4 „B 0-6 0-8 



y[rad] "m 

Figure 3. Marginalised Bianchi VIIj, parameters for left-handed model G 
using the original WMA P ILC map (red) and that connected for the texture 
fit o flCruzet alj j2007bl) (black). 



4 CONCLUSIONS 

In [Bridges et alj ( |2007|) we concluded that the preferred Bianchi 
Vllh template was likely to be heavily influenced by the cold spot 
in t he southern galact ic hemisphere. Accepting the recent texture 
fit ( Cruz et al.' 2007b) essentially removes this feature from the 
WMAP data. Although there were other discernible features that 
appeared reduced following 'correction' of the data by the best- 
fit Bianchi VIIi, template, these were mostly concentrated close 
to the galactic plane where little confidence can be placed in the 
map quality. The Bayesian evidence in favour of the inclusion of 
a Bianchi Vllh component, was only marginally significant for just 
two of the models we had previously considered. Using data 'cor- 
rected' for the texture, both of these models are now left almost 
entirely unconstrained by the data and consequently register dis- 
favouring log evidence values. This result raises further doubts over 
the relevance of a Bianchi Vllh explanation for any of the anoma- 
lous features seen in the current WMAP data. 
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